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Abstract. We present the results of a very large aspect ratio (A = 42.6) numerical simulation of fully compressible turbulent 
convection in a polytropic atmosphere, and focus on the properties of large-scale flows. Mesoscale patterns dominate the 
turbulent energy spec trum. We show that these structures, which had already been observed in Boussinesq simulations by 
ICattaneo et al. ( 2001), have a genuine convective origin and do not result directly from collective interactions of the smaller 
scales of the flow, even though their growth is strongly affected by nonlinear transfers. If this result is relevant to the solar 
photosphere, it suggests that the dominant convective mode below the Sun's surface may be at mesoscales. 
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1. Introduction 

The origin of solar photospheric flows on horizontal scales 
larger than granulation (C ~ 1 000 km) has been a puz- 
zling problem for more than forty years, when s upergr anu- 
lation (I ~ 30 000 km) was discovered by iHartl dl954h and 
later on confirmed bv lSimon & Leighton ( 1964b . Even though 
recent breakthroughs in the field of supergranulation imaging 
have been made thanks to the emergenc e of local helioseismol- 
ogy techniques dDuvall & Gizonl|2000h and the results of the 
MDI instrument dHathawav et alj|2000r) . its origin is still un- 
clear. The existence of an intermediate scale, mesogranulation 
(I ~ 8 000 km), is als o a matter of debate dHathawav et al 
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Meso and supergranulation have long been believed to be 
due to Helium deep recombinations driving cell- like convec - 
tion. This view now appears to be out of dat e dRast 20031) 



Several numerical experiment s of convection dCattaneo et al 



199 it IStein & Nordlundfe oOO) at moderate aspect ratio (A is 
the ratio of the box width to its depth) have shown a ten- 
dency of long-lived large-scale flows to form in depth. Using 
large asp ect ratio (A = 20) simulations of Boussinesq con- 
vection, ICattaneo et al.l d200lb have suggested that mesogran- 
ulati on may res ult from nonlinear interacti ons of granules (see 
also|Rast 2003). A large scale inst ability dGama et alJll994 ) of 
granules has also been proposed bv lRieutord et alj d2000l) to ex- 
plain supergranulatio n. Local numerical simulations at A = 1 
dRieutord et ai]|2002l) did not confirm it. beRosa et all d2002l) . 
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using spherical simulations, have computed flows down to su- 
pergranular scales. Actually, the emergence of the three distinct 
scales of granulation, mesogranulation, and supergranulation in 
the surface layers, among the observed continuum of scales, re- 
mains a fully open problem that still deserves much work. 

In this Letter, we report new results on three-dimensional 
numerical simulations of fully compressible turbulent convec- 
tion in a rectangular box with very large aspect ratio A = 42.6. 
This configuration allows us to study accurately the turbulent 
dynamics at horizontal scales between granulation and super- 
granulation, which have not been covered by previous numeri- 
cal simulations. A compressible fluid is used to provide a more 
realistic model of photospheric convection than a Boussinesq 
fluid. Also, density stratificati on should attenuate th e effect of 
an artificial bottom boundary (Nordlu nd et aljfl994 ). 

In Sect. |2] we present our numerical setup and physical 
model. Section [3] is devoted to the analysis of the flow. The 
main consequences of the results are discussed in Sect. @] 
which is followed by a short conclusion. 



2. Numerical model and run parameters 

For the purpose of our investigations we use a code designed 
to solve the fully c ompressible hydrod ynamic equations for a 
perfect gas (e. g. ICattaneo et al.ll 1 99 lb in cartesian geometry. 
Constant dynamical viscosity and thermal conductivity are as- 
sumed. A constant thermal flux is imposed at the bottom, while 
temperature is held fixed at the surface. The velocity field sat- 
isfies stress-free impenetrable boundary conditions. The ini- 
tial state is a m = 1 polytropic atmosphere (y = 5/3) with 
small random velocity v, temperature 6, and density p pertur- 
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bations. The initial density contrast between the bottom and 
top plates is 3, a value for which most of the features of strat- 
ification are already present in the linear convective instability 
problem (IGough et al.lll976l) . The Prandtl number is Pr = 0.3 
and the Rayleigh number evaluated in the middle of the layer 
is R = 3 1 5 (650 times supercritical). We do not paste sev- 
eral smaller boxes in itially, as was the case in the paper by 
Catt aneo et al.1 (l200lh . thus there is no artificial spatial sym- 
metry at t — 0. 

A sixth-order compact finite difference scheme ( Leldll992 ) 
is used in the vertical (gravity g) direction and a spectral 
scheme in the horizontal (periodic) directio ns. FFTs are im- 
pleme nted via the MPI version of FFTW (IFrigo & Johnson 
1998 |. Dealiasing b y removal is performed using the 2/3 rule 
( Canuto et al.ll 1 988b . Time-stepping is done with a third-order, 
low-storage, fully explicit Runge-Kutta scheme. Energy dissi- 
pation is handled by laplacian terms, without any subgrid-scale 
modelling or hyperviscosity. A very large aspect ratio A =42.6 
was achieved using 82 x 1024 x 1024 grid points. The simula- 
tion ran on 64 processors and a total of 400 GB of raw data was 
collected throughout the numerical experiment. 
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Fig. 1. Normalized momentum E(k, z) and temperature Eg(k, z) 
power spectra at t — 0.7 and z = 0.99 (the temperature spec- 
trum has been shifted). The maximum around k ~ 7 cor- 
responds to mesoscales. Also, a significant power excess in 
Eg(k,z) is observed around k ~ 40 close to the surface, in 
comparison to the deeper layers. This feature is associated with 
granulation in this simulation. 



3. Results 

In the following, the depth of the layer d is used as unit of 
length and the vertical thermal diffusion time d 2 /K\, ot as unit of 
time (/Cbot is the thermal diffusivity at the bottom of the layer). 

3.1. Flow structure and evolution 

We first describe the flow evolution during 0.7 thermal diffu- 
sion time, corresponding to twelve turnover times (twice the 
vertical crossing time based on the r.m.s. velocity). Initially, lin- 
ear growth is observed for a normalized horizontal wavenum- 
ber k ~ 45 (length I ~ 1) predicted by linear theory. The maxi- 
mum of the depth-dependent momentum spectrum (hats denote 
two-dimensional horizontal Fourier transforms) 

E(k,z)= f \{p) k (z)\ 2 kd£l k (1) 
Jn k 1 1 

then progressively shifts from k ~ 45 to k ~ 5 - 8 {( ~ 5 - 8) at 
the end of the simulation (Fig.[TJ. As can be seen on Fig. [2] the 
integral scale reaches L m ~ 7. Horizontal temperature maps at 
different times (Fig. [3} clearly show that the dominant scale of 
the flow increases during the simulation. This coherent pattern 
is best seen in the middle of the layer, but most of the in-depth 
dynamics are still clearly visible in the surface layers. In the up- 
per layers, a second distinct smaller scale appears (bottom-right 
picture of Fig . [3j . Its vertical extent (0.2 d) corresponds to the 
surface region of superadiabatic stratification, while the interior 
is almost isentropic. Referring to the Sun, we shall identify this 
thermal boundary layer scale with granulation, thus the larger 
internal scale is a mesoscale. 

Next we shall try to explain the growth and saturation of the 
integral scale as well as the origin of the observed mesoscales. 
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Fig. 2. Integral scale L int (z) = § J k~ l E(k, z) dk \ jE(k,z)dk 
evolution at z — 0.87. This scale represents the most energetic 
flow structures. 



3.2. Convective origin of mesoscales 

After the linear growth, saturation of the velocity amplitude 
at k ~ 45 occurs. At that time, this mode is the scale of 
energy injection. Larger scales are also highly linearly un- 
stable for R = 3 10 5 and we observe that they continue to 
grow. However, in the nonlinear regime, energy transfer to 
smaller scales limits this growth. To show this, we take the 
horizontal Fourier transform of the momentum equation and 
extract its solenoidal part by applying the projection opera- 
tor P[v] = v - VA~'V-v = v 5 , denoted by V when acting on 
horizontally Fourier- transformed fields. Taking the dot product 
with the complex conjugate of the solenoidal part of momen- 
tum (pv)\(z), integrating over depth z and angles in the 




Fig. 3. From left to right and top to bottom, time-evolution of temperature maps at z - 0.63. The size of the visible mesoscale 
pattern increases until a quasi-steady state is reached (bottom-left image). Bottom-right picture: same as bottom-left image, for 
a surface layer (z = 0.99), showing the differences between in-depth mesoscale dynamics and the smaller-scale flow in the upper 
thermal boundary layer. 



horizontal spectral plane, we obtain a time-evolution equation 
for the solenoidal part E s (k, z) of E(k, z), integrated vertically: 



Jo 



d, E\k, z) dz = T(k) + F(k) + D(k), 



T(k) = -2 Re 



F(k) = 2 Re 



f f (pv)Vp[V-(pvv)j£dn*dz 
Jo Jn t 

Jo Jn 



(pvy k -<p[pg k ]kdn k d z 



(2) 



(3) 



(4) 



T(k) represents nonlinear transfers, F (k) is the forcing by buoy- 
ancy and D(k), which has a similar definition, represents vis- 
cous dissipation. In Fig. HI we plot these quantities averaged 
over a time interval during which k > 12 modes are steady, 
while modes with smaller k develop. We observe that F{k) is 
the basic energy supply on large scales, as in the linear convec- 
tive instability mechanism. Nonlinear transfer T(k) is always 
negative for modes with 7 < k < 11, which have a small but 
positive net energy growth. Therefore large scales do not come 
out of nonlinear interactions but have a convective origin. This 
effect could also be observed using the energy equation: as in 
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Fig. 4. Depth-integrated spectral transfer T(k), buoyancy 
forcing F(k), dissipation D(k) and net power injected in 
jE s (k,z)6z following Eq. ©. An average between t = 0.14 
and t = 0.20 has been taken to outline the mean growth of 
k < 12 modes during this period owing to F(k) (note especially 
that T(k) < for these scales). The maximum over k of the 
depth-integrated spectrum has been used for normalizing. 
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linear convection, energy is injected in large scales via the ad- 
vection of the horizontally averaged entropy profile, while non- 
linear transfers and diffusion only remove energy from them. 

Note finally that the mesoscale pattern of Fig.[3]is expected 
to expand slightly on a much longer time scale that can not be 
achieved numerically. Also, the dominant scales may depend 
on the Rayleigh number, as they result from a balance between 
buoyancy and nonlinear transfers. The se mesocells are very 
probably the same as those observed bv lCattaneo et alj (1200 ll) 
in Boussinesq simulations with A = 20. Their size is compara- 
ble in both experiments. We therefore confirm these results for 
a compressible fluid, in a larger aspect ratio box with no initial 
symmetry, but interpret them quite differently. 



4. Discussion 

4.1. Relations with solar photospheric convection 

We now discuss the relevance of our results to the Sun. First, an 
estimate of the size (in km) of the typical structures of our sim- 
ulation can be computed to clarify the comparison with obser- 
vations. Identifying the thermal boundary layer thickness with 
the typical vertical extent of solar granulation (~ 150 km), we 
find our mesocells to be 6 000 km wide. 

This numerical experiment represents a highly idealized 
model of photospheric convection, even though it integrates 
density stratification. Differ ences with observatio ns or with 
more realistic simulations (IStein & Nordlundll 19981) are there- 
fore clearly expected. The most important one is the prominent 
peak at mesoscales in our power spectrum (Fig.[TJ, which might 
be due to boundary condition effects or to the absence of ra 



diative transfer in our simulations. As noted by ICattaneo et al 



(1200 ll) . granulation is directly related to the formation of a ther- 
mal boundary layer, so that changes in boundary conditions 
might have a strong impact on the contrast between mesogran- 
ular and granular flows. A l so, in radiative convection simula- 
tions bv lStein & Nordlu nd (2000) with an open bottom bound- 
ary, the dominant scale visually increases continuously with 
depth, which does not happen in our experiment. Besides, the 
absence of radiative transfer in our simulations makes it im- 
possible to defin e a r = 1 surface. Actually, the intensity 
map presented bv lRieutord et al. (2002) does not exhibit clear 
mesoscale intensity modulation whereas temperature maps at 
fixed depth do. Since the r = 1 surface does not correspond to 
a fixed depth, mesoscale convective flows may be present in the 
subsurface layers and be partly hidden this way. 

Besides , various ob servations suggest some mesoscale or- 
ganization. Odal ( 1984 ) has reported some clustering of gran- 
ules around brighter granules distributed on mesoscales. In 
our simulation (Fig. [3), the imprint of mesoscales is very 
clear at the surface and clustering of smaller granules also 
occurs around bright spots corresponding to mesoscale up- 
flows. Finally, a distribution of inter-network magnetic fields 
at the scale of mesogranulation has been found recently 



4.2. Very large scales 

A final word must be said about very large-scale flows. 
According to the previous estimates, the horizontal size of our 
box is 35 000 km, which does not leave much room for su- 
pergranules. We do not obtain a second peak at small k in the 
spectra, so there is strictly speaking no supergranulation in our 
simulation. This may be because we do not have the neces- 
sary physical ingredients in the model, because our box is not 
large or deep enough, or because our run is too short. However, 
noticeable positive nonlinear transfer occurs at the smallest k 
(Fig. H]i and we observed weak horizontally divergent large- 
scale flows of still unclear origin. When subjected to horizontal 
strains, strong mesoscale vertical vortices resulti ng from angu- 
lar m omentum conservation in sinking plumes (IToomre et alJ 
1990) might help creating horizontal vorticity on large scales. 



(DormnguezCerdena 2003; Roudier & Muller 2004), which 



might also be related to strong convective mesoscale plumes. 



5. Conclusions 

The results of our simulations of large-scale convection in a 
compressible polytropic atmosphere show that the dominant 
convective mode is found at mesoscales. We observe that large 
aspect ratio simulations are necessary to study the convective 
dynamics of these structures, since the integral scale is L; nt ~ 7 
in the final quasi-steady state. A slow evolution is expected on 
longer time scales, which are unfortunately numerically out of 
reach. Some similarities with solar observations are found. If 
this kind of model is relevant to study the solar photosphere, 
our results suggest that mesoscale convection may be powerful 
below the Sun's surface. This would help to explain the coex- 
istence of two apparently distinct granular and mesogranular 
scales. Supergranulation is not found in this experiment. 
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